Coulomb corrections in the extraction of the proton radius 
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Multi-photon exchange contributions are important in extracting the proton charge radius from 
elastic electron-proton scattering. So far, only diagrams associated with the exchange of a second 
photon have been evaluated. At the very low Q 2 values relevant to the radius extraction, higher order 
contributions may become important. We evaluate these corrections in the Effective Momentum 
Approximation, which includes the Coulomb interaction to all orders, and find small corrections 
with a strong Q 2 dependence at low Q 2 and large scattering angles. This suggests that the higher 
order terms may be important in the evaluation of the proton magnetic radius. 

PACS numbers: 13.40.Gp,14.20.Dh,25.30.Bf 



The proton RMS charge radius, R p , has become a topic 
of great interest following the recent Lamb shift measure- 
ment in muonic hydrogen [l[ which is extremely sensi- 
tive to R p . This measurement yields R p —0. 8418(7) fm, 
significantly smaller than recent extractions based on 
electron-proton interactions 043- These include both 
electron scattering and electron-proton interactions in 
the hydrogen atom, which yield a combined result of 
i? p =0.8772(46) [!, Q. Replacing the CODATA06 value 
from Ref. H with the updated CODATA10 value 
yields an electron-based average of 0.8775(39), corre- 
sponding to a 9cr difference between the muonic hydrogen 
extraction and the combined electron measurements. 

Radiative corrections play an important role in the 
electron scattering measurements. While the largest cor- 
rections are well understood, the diagrams which depend 
on the proton structure, e.g. the two-photon exchange 
(TPE) diagrams, cannot be calculated exactly due to 
the hadronic structure uncertainty. In the past, these 
corrections were calculated in the 2nd Born approxima- 
tion, assuming the exchange of a second soft photon with 
an unexcited intermediate state. Initial calculations were 
performed in the limit Q 2 — > Q, and later at finite Q 2 , 
which require a parameterization of the proton charge 
and magnetic form factors, Ge and Gm- The inclusion 
of these corrections was found to be important in the 
extraction of the proton radius 0, [j| . 

More recently, calculations going beyond the 2nd Born 
approximation have been performed [lOMlq . motivated 
by the discrepancy between Rosenbluth and polariza- 
tion measurements at high Q 2 17H2ll|. See recent re- 
views 22J, |23j for further details on the different theoreti- 



cal approaches. There have also been seve ral p henomeno- 
logical extractions of TPE contributions |2J-|30j and at- 
tempts to experimentally constrain TPE j3ll - l37| . but es- 



sentially all of these works focus on Q > 2-3 GeV", 
where there is a clear discrepancy between the Rosen- 
bluth and polarization extractions of the form factors. 

Focusing on the low Q 2 regime, Figure [T] shows the 
range of TPE corrections to the cross section based on a 
set of calculations which include both hard and soft TPE 
and which aim to be reliable at low Q 2 : The hadronic 
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FIG. 1: (color online) R ang e of results for several TPE cal- 
culations at low Q 2 [n], Gj, 03, |H, [H. Solid, long-dashed, 
short-dashed, and dotted lines show the range for Q 2 =0.01, 
0.03, 0.1, and 0.2 GeV 2 , respectively. Note that the low Q 2 
expansion [38|] is only valid up to Q 2 =0.1 GeV 2 , and so is 
excluded from the Q 2 = 0.2 GeV 2 ranee. 



calculation of Blunden, Melnitchouk, and Tjon [111], a 



similar calculation by Borisyuk and Kobushkin 13j , their 
low-Q 2 TPE expansion 38], and their dispersion calcula- 
tion excluding [3] and including [3^| A resonance contri- 
butions. While the correction decreases with increasing 
Q 2 for e > 0.4, the TPE correction for e < 0.3 first in- 
creases and then decreases. These corrections have been 
shown to have an important effect on the form factor 
extraction at low Q 2 values [13], as well as impacting 
the extracted charge and magnetization radii of the pro- 
ton 0, [l3|, 41 - 43 1 . While the calculations use different 
approaches and have different input, they are in good 
agreement with a typical spread between calculations be- 
low 0.05%, suggesting that the TPE corrections are well 
known in this region. At these low Q 2 values, deviations 
from a linear correction in e are small except for extreme 
e values, making tests of the linearity of the Rosenbluth 
separation [3, [32| fairly insensitive to the low Q 2 contri- 
butions of these calculations. 
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Extraction of the radius requires cross sections mea- 
sured at low Q 2 , meaning low electron energies, espe- 
cially in the case of large-angle measurements needed 
to constrain Gm and the magnetization radius. At low 
energies, the change in the electron energy due to its 
interaction with the Coulomb potential of the nucleus 
may yield an important change in the kinematics at the 
scattering vertex. In quasi-elastic scattering this effect 
is sometimes estimated in the Effective Momentum Ap- 
proximation (EM A) where the acceleration due to 
the Coulomb field of a high-Z nucleus can have a sig- 
nificant impact on the e-p scattering kinematics. This 
approach has the advantage of including the Coulomb 
interaction to all orders in a semi-classical picture, and 
can easily be applied to both high- and low-energy kine- 
matics. This is to be compared to the 2nd Born approx- 
imation which is a more rigorous approach, but which 
includes only two-photon exchange. 



change observed when the 2nd Born calculation is applied 
to the data [§]. 

In this paper, we estimate the impact of higher-order 
Coulomb corrections by applying the EMA prescription 
of Ref. [45[ to elastic e-p scattering at low Q 2 . The key 
parameter in the calculation is the Coulomb potential at 
the point where the scattering occurs. When the EMA 
is used to evaluate scattering from a heavy nucleus, it is 
generally assumed that the scattering occurs uniformly 
within the nucleus, so the Coulomb potential is taken to 
be something between the surface and central potential. 
Averaging over the nuclear volume, assuming a uniform 
charge sphere, yields a potential that is 80% of the central 



potential, in good agreement with the result of Ref. [44 1 
which adjusts the average potential in the EMA calcu- 
lations to match distorted-wave Born approximation re- 
sults 0, |4?J which are more reliable but can only be 
performed at lower Q 2 values. 
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FIG. 2: (color online) Coulomb corrections calculated in the 
2nd Born approximation for e-p elastic scattering for Q 2 val- 
ues from 0.01 GeV 2 (top curve) to 1 GeV 2 (bottom curve). 

Figure [2] shows the fractional correction to the cross 
section in the 2nd Born approximation. One can make 
a crude estimate of the impact of these corrections on 
the radius. The RMS radius, R p , is defined in terms 
of the low Q 2 expansion of the form factor, Ge(Q 2 ) ~ 
1 — Q 2 R 2 /6, with a similar expression for the mag- 
netic radius. A proton RMS radius of 0.85 fm yields 
G E {Q 2 ) « 1 - 3Q 2 , with Q 2 in GeV 2 , or a fractional 
slope at Q 2 = of roughly 300%/GeV 2 . The charge 
form factor contribution to the cross section goes like 
G 2 E (Q 2 ), for a slope of roughly 600%/GeV 2 in the cross 
section. Coulomb corrections in the 2nd Born approxi- 
mation change the cross section by about 0.2% between 
Q 2 =0.01 and 0.03 GeV 2 (Fig. [2]), yielding a change in the 
slope of 0.2%/0.02 GeV 2 =10%/GeV 2 . This is roughly 
2% of the total slope introduced by the protons size, cor- 
responding to a 1% change in the extracted radius, in 
good agreement with the observation of a roughly 0.01 fm 



FIG. 3: (color online) Coulomb corrections calculated in 
the Effective Momentum Approximation (as described in the 
text) for e-p elastic scattering for 0.01 < Q 2 < 1 GeV 2 . 

At high energies, the EMA calculation yields a correc- 
tion that is very similar to the 2nd Born approximation 
for e > 0.5, as shown in Ref. 45[, with differences the 1- 
2% level for larger-angle scattering. Note that the EMA 
calculation applied in Ref. 45J used the central potential 
and did not include the flux factor of Ref. 44J . Inclusion 
of this correction yields somewhat improved agreement 
at low e, but there are still clear differences remaining. 

For e-p elastic scattering at very low Q 2 , one expects 
that the scattering will occur when the electron is outside 
of the proton, so using the average Coulomb potential 
in the proton would be an overestimate. We take the 
Coulomb potential corresponding to the case where the 
scattering occurs at a separation Ax = l/q, where q is 
the momentum of the exchanged virtual photon. This 
will suppress the impact of the energy shift for Q 2 values 
where the scattering occurs outside of the proton. This 



3 



occurs for Q 2 < 0.03 GeV 2 if we take the proton to be a 
uniform sphere of radius 1.15 fm to match the observed 
RMS radius. For Q 2 values where the scattering occurs 
inside of the proton, we limit the Coulomb potential to 
a maximum of the volume-averaged value of 1.5 McV. 

Figure [3] shows the Coulomb correction in the EMA. 
The behavior is qualitatively similar to the 2nd Born ap- 
proximation result for Q 2 — 1 GeV 2 , as it is for larger Q 2 
values 451 ]. The EMA result has a very different angu- 
lar dependence at very low Q 2 values, with a sharp rise 
at low e. As such, this correction could have an impor- 
tant impact on the low Q 2 extraction of the form factors 
and radius, especially for the extraction of Gm and the 
magnetic radius. 
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FIG. 4: (color online) Q 2 dependence of the EMA correction 
at £=0.02 (solid), 0.2 (long dash), and 0.5 (short dash). 

Figure 0] shows the EMA result as a function of Q 2 
for three e values. At larger e values, the effect is small, 
especially near Q 2 =0. For e = 0.02, corresponding to 
scattering near 180 degrees, there is a very rapid rise 
with Q 2 (note that Q 2 is shown on a logarithmic scale). 
From the Q 2 dependence of the correction, it is straight- 
forward to estimate what impact this would have on a 
direct extraction of the magnetic radius from very low 
Q 2 data at e w 0. The EMA yields a 1% change in the 
correction for a change in Q 2 of less than 0.01 GeV 2 , or a 
change in the slope of 100%/Q 2 . This would correspond 
to a very large correction to the magnetic radius. 

However, existing measurements below Q 2 =0.1 GeV 2 
are limited and there are not measurements over this Q 2 
and e range. The kinematics of interest are at some- 
what larger Q 2 or e values. At £=0.02, the correction de- 
creases by 1.5% going from Q 2 = 0.03 GeV 2 to 0.2 GeV 2 , 
while for e = 0.2, it increases 0.3% in going from 0.01 to 
0.03 GeV 2 . These provide changes in the fractional slope 
of the cross section of -8%/GeV 2 and +15%/GeV 2 , re- 
spectively, corresponding to small but potentially signifi- 
cant corrections to the extracted magnetic proton radius. 
The corrections at larger Q 2 and e values are smaller, but 



these kinematics require extrapolations to e = and low 
Q 2 , and thus even these smaller corrections may not be 
entirely negligible. 

While this rough estimate shows that these effects 
could have important contributions to the extraction of 
the magnetic radius, it is difficult to provide a more de- 
tailed estimate of the impact. First, because it is very 
sensitive to the exact kinematics of the data included in 
the extraction. The corrections would have to be applied 
to each measurement and then the radius extraction pro- 
cedure, including fitting of normalization factors if they 
are allowed to vary, would have to be repeated. In ad- 
dition, it is not clear how one would combine the results 
of the EMA at very low Q 2 with the TPE calculations, 
to include the full impact of the corrections without dou- 
ble counting. A more detailed calculation, e.g. in the 
distorted-wave Born approximation, would allow for a 
more reliable evaluation of the importance of the higher 
order contributions. 

In conclusion, we find a clear difference between the 
2nd Born approximation and the EMA evaluation of 
Coulomb corrections for e-p elastic scattering at low 
Q 2 . The difference is particularly important large an- 
gles, where the data are sensitive to Gm(Q 2 ), suggesting 
that effects beyond two-photon exchange may be impor- 
tant in the extraction of the magnetic radius. A realistic 
estimate of their impact would involve a more detailed 
calculation, with the impact of the correction evaluated 
using on the exact data set and fitting procedure used to 
extract the radius. 
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